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Abstract. We show that dark matter abundance and the inflationary scale H could be 
intimately related. Standard Model extensions with Higgs mediated couplings to new physics 
typically contain extra scalars displaced from vacuum during inflation. If their coupling to 
Standard Model is weak, they will not thermalize and may easily constitute too much dark 
matter reminiscent to the moduli problem. As an example we consider Standard Model 
extended by a Zo symmetric singlet s coupled to the Standard Model Higgs <f> via A4A<f>,s 2 . 
Dark matter relic density is generated non-thermally for A < 1CU'. We show that the dark 
matter yield crucially depends on the inflationary scale. For H 10 10 GeV we find that the 
singlet self-coupling and mass should lie in the regime A s > 1CU 9 and m s < 50 GeV to avoid 
dark matter overproduction. 

Keywords: Dark matter, Inflation, Freeze-in, SM Higgs, Higgs portal 


Contents 


1 Introduction 1 

2 Singlet dark matter and freeze-in without primordial condensates 2 

3 Primordial condensates 3 

3.1 Inflationary fluctuations and their evolution 3 

3.2 The moduli problem 4 

4 Freeze-in with the primordial condensates 5 

4.1 Processes in the quartic regime 6 

4.2 Processes in the quadratic regime 7 

5 Estimating the total dark matter yield 7 

5.1 Case I: The condensate decays completely 8 

5.2 Case II: The condensate survives 9 

6 Conclusions and outlook 11 


1 Introduction 

For the measured Higgs mass — 125 GeV the Standard Model (SM) vacuum is metastable 

up to remarkably high energies /r c ~ 10 1 * * * * * * * * * 11 * * GeV with a lifetime much longer than the age of 

the universe [1-4]. In the early universe the stability is crucially affected also by the large 

spacetime curvature. Radiative corrections to the Higgs potential necessarily induce the 

non-minimal curvature coupling which during inflation when R ~ H 2 may easily 

dominate the potential [5-7]. Indeed, it has been shown that if £(/tew) >0.1 [5] at the 

electroweak scale the SM Higgs remains stable up to the highest inflationary scale 77* ~ 10 14 

GeV consistent with the tensor bound [8, 9] without any need to include new physics. 

New physics beyond SM is however implied by several cosmological observations, such 

as inflation itself, dark matter and baryon asymmetry. If the Higgs potential is not drasti¬ 

cally modified by new physics the Higgs generically is a light and energetically subdominant 
field during inflation [5, 6, 10-17]. (The attracting possibility of Higgs-driven inflation re¬ 
quires a non-trivial deviation from the SM potential at high energies, such as a large coupling 
to spacetime curvature [18-23].) Inflationary fluctuations displace the light Higgs from its 
vacuum generating a primordial Higgs condensate [10, 11]. The resulting out-of-equilibrium 
initial conditions /i* ~ 77* for the hot big bang epoch may have significant observational 
ramifications. Particle production from the time-dependent Higgs condensate may for ex¬ 

ample generate baryon asymmetry [24] or produce non-thermal dark matter [25]. A careful 
investigation of the observational effects is of key importance and could reveal powerful new 
tests of specific SM extensions at very high energy scales. 

In this work we will investigate how the initial conditions set by inflation affect the 
generation of dark matter abundance in the interesting class of portal scenarios [26-28] where 
the SM fields feel new physics only through Higgs-mediated couplings. As a representative 
example we will consider a Z 2 symmetric scalar singlet s coupled to the Higgs field by 
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Vint = A 2 h 4d<I>s 2 . This simplified example captures some interesting features of the portal 
models. The singlet constitutes a dark matter candidate [26, 27, 29, 30] and its dynamics 
at electroweak phase transition could allow for baryogenesis [31-33]. Analogously to Higgs 
also the singlet is generically a light field during inflation and gets displaced from vacuum 
[25]. In this simple setup the out-of-equilibrium initial conditions do not typically affect the 
electroweak baryogenesis as the scalar condensates have either decayed or diluted away by 
that time [25]. They could however significantly affect the dark matter abundance if the 
portal coupling is very weak A s h < 1CU 7 . In this case the singlet never thermalizes with 
the SM fields. The singlet dark matter originates entirely from non-thermal production of 
singlet particles through the so called freeze-in mechanism [27, 34-43]. As we will show in this 
paper, the presence of scalar condensates significantly alters the previous estimates for the 
efficiency of the process leading to novel interplay between dark matter properties and initial 
conditions sensitive to the inflationary scale. Moreover, the inflationary displacement of the 
singlet may also lead to overproduction of dark matter which places stringent constraints on 
viable singlet mass scales and the values of its self-coupling. 

The paper is organized as follows: In Section 2 we define the Z 2 symmetric model and 
discuss the freeze-in mechanism of dark matter production. In Section 3 we discuss initial 
conditions set by inflation and the related moduli problem. In Section 4 we investigate 
the production of singlet dark matter through the freeze-in accounting for the presence of 
primordial scalar condensates in the Boltzmann equations. We then present our main results 
for the dark matter abundance and its dependence on scalar couplings and the inflationary 
scale in Section 5. Finally we summarize and discuss the results in Section 6. 

2 Singlet dark matter and freeze-in without primordial condensates 

We consider the SM extended to include a Z 2 symmetric scalar singlet [26, 29, 30] 

U(4>, S ) = + Ah(4» t 4>) 2 + ^m s V + js 4 + ^($ t <h)s 2 . (2.1) 

Here 4) is the Standard Model Higgs doublet and s is a real scalar singlet assumed to possess 

a Z 2 symmetry in order to make it stable. Therefore, the singlet constitutes a dark matter 

T 

candidate. At low temperatures and in the unitary gauge, 4> = (0, (u + h)/y/ 2) . We will 
assume that to 2 > 0 and A s h > 0, so there is no spontaneous symmetry breaking in the singlet 
sector. This guarantees that the singlet scalar is stable in the vacuum. Moreover, in the cases 
we will consider, the portal coupling A s h is assumed to be very weak. This guarantees that 
for light enough singlets there are no constraints from the invisible decay width of the Higgs 
at the LHC. 

Usually the exact value of the self-interaction A s is considered to be irrelevant for dark 
matter abundance. However, as we shall see it is of uttermost importance in determining the 
initial conditions for low energy phenomena and, consequently, for the calculation of the total 
dark matter yield via the freeze-in mechanism relevant in the limit of weak portal coupling 

Ash < nr 7 . ” 

To set the notation, we start with a brief review of the well-known freeze-in results 
[27, 34-43] obtained neglecting the impacts of primordial scalar condensates generated by 
inflation. 

In the freeze-in setup the singlet dark matter is produced from the thermal bath of 
SM particles through out-of-equilibrium decays and scatterings. For concreteness, consider 
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the case where the dominant process is Higgs decay into two singlet particles [34] which is 
possible below Tew whenever > 2m s . For a discussion of other Higgs mediated processes 
which for < 2 m s could also be important, see e.g. [35]. 

The evolution of number density of the singlet scalar is determined by the Boltzmann 
equation 


h s + 3Hn s = J dU h dU Sl dU S2 (2TT) 4 6 4 {p h - p Sl - p S2 ) 

X {\M\U ss f h (l + / s )(l + f s ) - \M\ 2 ss ^ h f s f s (l + / h )) , (2.2) 

where cTIi = d 3 ki/((2ir) 3 2Ei), M. is the transition amplitude and /; is the usual phase space 
density of particle i. The Higgs particles are assumed to be in thermal equilibrium, and in 
the usual approximation one assumes that Maxwell-Boltzmann statistics can be used instead 
of Bose-Einstein, ~ e ^ Eh l T . 

Setting f s = 0 on the right hand side of Eq. (2.2) the singlet abundance, produced at 
low temperatures by thermal Higgs particles only, then becomes [34] 


P s /i 2 « 1.73 x 10" 


27 m s^h->ss 

9 

mt 


= 1.73 x 10 27 ^ 


m s 


■ (2 ' 3) 


mj] \?>2'Km\ l 


In the limit, m s <C m h, this yields a parametric estimate for the coupling sufficient to produce 
a sizeable dark matter abundance 


A sh - 10 


-li 


'n s h 2 \ 1/2 f GeV A 1/2 


mrc- 


(2.4) 


The implied small coupling values are compatible with the key assumption of the freeze-in 
scenario that the dark matter candidate does not thermalize with the SM background above 
the EW scale. 

The essential approximation in the above analysis, and commonly made in all freeze-in 
computations, is that the dark matter abundance is initially negligible, i.e. f s = 0 on the 
right hand side of Eq. (2.2) [34, 35, 38]. However, this may not be the generic outcome 
in realistic setups where ramifications of the infiationary stage are consistently accounted 
for. Indeed, scalar dark matter candidates could easily get displaced from vacuum during 
inflation and their subsequent relaxation towards the low energy vacuum can significantly 
alter the freeze-in picture as we will now turn to discuss. 


3 Primordial condensates 

3.1 Inflationary fluctuations and their evolution 

During inflation any energetically subdonrinant light scalar acquires superhorizon fluctua¬ 
tions proportional to the inflationary scale 77*. The spectators get locally displaced from 
their vacuum state and after the end of inflation the observable patch generically features 
primordial spectator condensates. 

Within the extended SM (2.1) we are considering here, both the Higgs and the singlet 
are light spectators during inflation, V" <C H 2 , and get displaced from vacuum [11, 25]. The 
typical amplitudes of the primordial Higgs and singlet condensates are given by the root 
mean square of fluctuations which for A s h < \/A s Ah yields [44] 

K=0( 0.1)^, s, = 0( 0.1)^. (3.1) 
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The inflationary scale H * can be expressed in terms of the tensor to scalar ratio r and 
amplitude of the scalar perturbations V £ as 


27T 7 


= -Vc 


(3.2) 


In the following we adopt these generic order of magnitude estimates as the initial conditions 
for the hot big bang epoch and explore their impacts on the singlet dark matter yield. 

Here we have assumed that the inflationary scale is below the flat space instability scale 
77* < 10 11 GeV such that A(//) > 0 and the curvature induced effective Higgs mass £Rh? is 
negligible. If the inflationary scale is higher, one should carefully account for the non-minimal 
curvature coupling which dominates the Higgs dynamics when A —> 0 [5, 6]. 

Assuming an efficient reheating of the SM sector soon after the end of inflation the Higgs 
condensate will be rapidly destroyed by the thermal bath [25] (see [11, 45, 46] for the decay 
at zero temperature and [47, 48] for reheating in the context of Higgs inflation). The singlet 
condensate on the other hand will not feel the thermal bath for A s h < 10~ 7 . Its evolution is 
affected by the redshifting due to expansion of space and also by out-of-equilibrium decays 
into singlet particles and Higgses through the portal coupling. 

Neglecting the decay processes for the time being, the condensate stays nearly constant 
until H 2 ~ A s s* after which it starts to oscillate with a decreasing envelope 


s 0 {T) ~ 


10 -3 A s 3/ V/ 4 t , 

10- 4 A s " 5/16 r 3 / 8 m s - 1/2 T 3 / 2 , 


T > T trans = 200A S 1/8 r -1 / 4 m s 

T < T 
L r^j ^trans • 


(3.3) 


For temperatures above T trans the singlet sees an effectively quartic potential A s s 4 m s s 2 
and its energy density scales as radiation, p s oc a~ 4 . Below T trans the quadratic mass term 
takes over and the singlet energy density scales as non-relativistic matter, p s oc a -3 . 


3.2 The moduli problem 

If the weakly coupled singlet condensate enters the regime p s ~ a -3 before the matter- 
radiation equality its energy density may lead to overproduction of cold dark matter. This is 
essentially the well-known moduli problem encountered in a variety of different theories for 
early universe physics (see e.g. [49-51]). 

In the SM extension by a scalar singlet (2.1) the moduli problem constrains the viable 
singlet mass scale m s from above. Using (3.3) and requiring that the energy density of the 
singlet condensate at matter-radiation transition T eq ~ 0.8 eV does not exceed the dark 
matter contribution I2 S < IIdm we obtain the mass bound 


m s 

GeV 


< 10 


-5 


A s 

io - 10 



(3.4) 


This constraint should be regarded as an absolute upper bound as it neglects the decay 
processes of condensate. The decays alleviate the constraint by causing the condensate 
amplitude to decrease faster and by depleting the singlet energy into the SM sector through 
Higgs mediated processes. 
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4 Freeze-in with the primordial condensates 


The primordial singlet condensate will significantly alter the standard freeze-in picture [34, 
35, 38]. Due to the time dependent background, singlet particles can be produced already 
well above the electroweak scale where h = 0. In this regime the effective Boltzmann equation 
for the number density of singlet particles takes the form 

h s + 3Hn s = j dPS S0tSMM2 \M\ 2 h ^ sh f s J h (l + f s )(l + h) 

+ j dPS S0jSl ' S2jS3 \M\^ssfs 0 fs(l + /s)(l + /s) 

/so/s(l + /h)(l + /h) 

+ r§Q— yss'W'SQ • (4.1) 


Here dPS SQ ^,c denotes phase space measure for the condensate so and particles a,b,c and 
contains the four-momentum conserving delta function. Furthermore n SQ denotes the con¬ 
densate number density, 


n S0 


d 3 k = p S0 

(2tt) 3 Js ° “ m Sj eff 


(4.2) 


where the singlet effective mass is defined as Tn'j, = V". As the coherently oscillating back¬ 
ground can only decay, inverse processes where particles would go back into the condensate 
are not present. 

The amplitudes M a _>.ft c in (4.1) correspond to decay processes induced by the oscillating 
condensate and denotes the rate for particle production directly from the time 

dependent effective potential. We approximate the amplitudes by [52-54] 


2ir5 4 (p 2 - pi)M i ^2 = 


/ OO 

dt(2\V(t)\l) , 

-oo 


(4.3) 


where V(t) denotes the interaction Hamiltonian induced by the singlet condensate 1 . For 
example, in the case of s h + h we have 


V(t) 


—A s hSo(^) / d?xshh 


(4.4) 


To extract the leading contribution for the particle production induced by the conden¬ 
sate, we linearize the Boltzmann equation in /}, and f s . After this the phase space integrals 
can be performed and (4.1) reduces to 


h s + 3 Hn s ~ 


K X (^) 

I<2 (^) 


F h^fsh n \\ T r s 0 — 


So ? 


(4.5) 


where K n is the nth modified Bessel function of the 2nd kind. On the right hand side we 
have neglected the source terms related to the processes s —> ss and s —>• hh. These are 

1 In the quartic regime, A a s 4 m B s 2 , the system is conformal and the amplitude coincides with the 

Minkowski result. In the quadratic regime we keep using the Minkowski metric neglecting the small 0(H/m s ) 

curvature corrections during one oscillation cycle. 
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suppressed by the small occupation numbers of singlet particles, / s <C /h, / So , and the latter 
are also kinematically heavily suppressed in the quartic regime. 

The Boltzmann equation for particles (4.1) is accompanied by the corresponding equa¬ 
tion of motion for the singlet condensate 


n So + 3 Hn So 




(4.6) 


The three contributions on the right hand side correspond to energy loss due to production 
of singlet particles and Higgses out of the oscillating condensate. The decay processes have 
negligible effects on the condensate motion until T ~ H. Up to this point the background 
dynamics is therefore well described by the solution (3.3). As T ~ H the amplitude of 
the condensate starts to decrease exponentially and to sufficient accuracy we can model the 
process as an instant decay at T = H. Consequently, the source terms in (4.5) vanish and the 
comoving singlet number density a 3 n s freezes to a constant. At the electroweak transition 
the generation of Higgs vacuum expectation value (vev) induces an additional contribution 
to the singlet number through the standard freeze-in mechanism. 


4.1 Processes in the quartic regime 

The decay rates in the Boltzmann equation take different values depending on whether the 
singlet is oscillating in the quartic So <1 rn s/ V^s or quadratic So < m s /yf\ ~ s regime of its 
potential. In terms of the temperature, these regimes correspond respectively to T > T trans 
and T < T trans , where the transition temperature is given by Eq. (3.3). 

In the quartic regime T > T trans , the amplitude of the oscillating condensate scales 
as so oc T. The rate for the h —> sh process in Eqs. (4.5) and (4.6) computed in the 
time-dependent background is given by 


p(a) _ 41 y- 5 tl 2 % 

*■“* 24lr h -Bf + El, 


= io- 4 a s 1/2 a1|| . 


(4.7) 


Here ky is the final state momentum satisfying the energy conservation condition E\ q + = 

mh + nw, where ui ~ 0.489m s ^ e ff is the oscillation frequency of the singlet condensate [55], 
and 

OO 

so(t)= s n& +iunt ■ (4.8) 

n =—oo 

The other processes in Eqs. (4.5) and (4.6) correspond to transitions from vacuum 
to two singlet or two Higgs states, induced respectively by the interactions A s SQ(f)s 2 and 
Ash'So(i)^ 2 - The corresponding rates are given by 


rW „ -4 X 10- 2 ^ S0 


r 


(4) 


A 2 


sh m ^eS 


S0 ^ hh 16yr 
where ?n 2 s = 3A s (so), and 


Ps 




'1- 


n =1 


2mh 

n<f> 


~ 0, 


So(*) - ( S 0> = Cn 


0 —i(f)nt 


(4.9) 

(4.10) 


(4.11) 
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with cj) being the oscillation frequency of £(f). The rate (4.10) is negligible due to kinematical 
suppression by the thermal Higgs mass mh(T) 3 > 


4.2 Processes in the quadratic regime 

In the quadratic regime the production of singlet particles directly from the time dependent 

( 2 ) 

background is energetically forbidden and consequently = 0. The transition h —> sh 

induced by the time-dependent condensate is also kinematically blocked in this regime. The 
only possible process in Eqs. (4.5) and (4.6) is then the production of Higgses out of the 
time dependent background sq —> hh. The rate for this process is given by 


r (2) = 3lAJi ( mh Y 

s 0 ^>-hh 54 ^ mg y \m s J 


(4.12) 


We reiterate that we are considering transitions in the time dependent background which 
amounts to the kinematical condition different for example from the standard 1 —> 2 decay 
in vacuum. 

From Eqs. (3.3) and (4.12) we obtain in the radiation dominated background the result 

p(2) p(2) 

So —yhh so —^hh 

H ~ H 

where ftrans denotes the time of transition from quartic to quadratic oscillations. Therefore, 
we immediately see that an eventual decay of the condensate needs to take place before the 
onset of quadratic oscillations. If the condensate has not decayed by ftrans the subsequent 
decay rates remain negligible and the condensate survives undecayed. 


Q-trans 

a 


(4.13) 


5 Estimating the total dark matter yield 

The primordial singlet condensate crucially alters the standard freeze-in estimates for the 
abundance and properties of singlet dark matter. Depending on the strength of singlet 
self-coupling, the condensate may either completely decay into singlet particles or survive 
comprising a coherently oscillating dark matter component. The two cases could have differ¬ 
ent ramifications on structure formation, see e.g. [56-59] . The fate of the singlet condensate 
also affects the eventual decay channels of the singlet provided that the portal sector contains 
additional degrees of freedom such as fermions with Yukawa couplings £ = gsi/jv/i. 

In addition to the primordial condensates, singlet particles are produced at the elec- 
troweak transition through Higgs decays A s h vhss. However, as opposed to the standard 
freeze-in estimates the singlet occupation numbers need not be small which could signifi¬ 
cantly affect the process. We leave a detailed analysis on this question to future work. In 
what follows we concentrate only on the abundace of singlet dark matter generated through 
the primordial condensates. 

We also note that the observational bounds are significantly different depending on 
whether the singlet constitutes isocurvature or adiabatic dark matter. While the component 
sourced by the primordial condensates clearly is isocurvature, the situation is less clear when 
the production of singlet particles through Higgs decay is important. Any additional cou¬ 
plings between the SM and the portal sector would also affect the situation. We will not 
dwell further on this important issue but simply choose to present both the adiabatic and 
isocurvature bounds in what follows. 
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5.1 Case I: The condensate decays completely 


As we have seen above, the condensate can decay only if the decay takes place in the quartic 
regime where A s Sq 3> m 2 . With the rates given in Eqs. (4.7) , (4.9) and (4.10) the Boltzmann 
equations for the particles (4.5) and condensate (4.6) in this regime are given by 


K x (^) 

n s + 3Hn s = — lmi , T h ^ sh n h + 

•^Sq— yss'fl'SQ ? 


h SQ + 3 Hn Sn = - 


so 


k 2 {W 

K\ (^) 


^ h-^fshJl h r 


so - 


(5.1) 

(5.2) 


Using that n SQ = p s /^s,eff = A s Sq/(4\/3) we can express the so dependence of 
(4.7) in terms of n So and rewrite the corresponding term in the condensate Boltzmann equa¬ 
tion (5.2) as 


„ _4*1 (=?)„. ft, 2 »h_ _f, „ 

^ mi, ^ ^h->-sh n h — 10 ^ 4v 3A sh ^ n2 n S0 _ Th->- s hn so . 

The solution for equation of motion (5.2) can then be written in the implicit form 


(5.3) 


^OSC 

a 


exp 



Sq^-SS Oo) ^h—>sh 


H 


H / 


(5.4) 


where the quantities with subscript osc are evaluated at t osc denoting the onset of quartic 
oscillations, H osc ~ A s s 2 and only the first term inside the integral depends on n so . 

When r so _>. ss , Fh-tsh < H the particle production clearly has a negligible effect on 
the motion of the condensate. Using that 2 r so _> ss oc s 0 , and r h-tsh oc T and substituting 
so oc a -1 under the integral in Eq. (5.4) the integral can then be performed. This yields the 
approximative explicit solution 


„ ( aosc Y ( r ^o^( g o) 

° _ A s 1/4 ' a ' ( H H 

Here i7* denotes the inflationary Hubble rate at the horizon exit of observable modes. 

For the weak portal couplings A s h < 1CU 7 we are considering here, we generically have 
r so ^ ss Th^ s h if the condensate is to decay in the quartic regime. The decay time is then 
dictated by r so ^ ss and the condensate decays practically instantaneously at i^dec = rs 0 -»-ss 
depleting all its energy into singlet particles. The solution of Boltzmann equation for singlet 
particles then reads 



n. 


1 ( K x (gjj T h ^ sh 

a 3 l k 2 (^) r S0 _> M 


n^a + n so a' 


t=t a 


(5.6) 


showing that the comoving number density freezes to a constant value as the time-dependent 
background vanishes. The behaviour is illustrated in Figure 1. 

The distribution of the generated singlet particles is peaked at k^/a^ec ~ V^s s o- As 
\/3A s so > m s , after the decay of the condensate the produced singlet particles constitute 

2 See Eqs. (4.9), (5.3) and (3.3), respectively, for the following scalings. 


- 8 - 













In a ’ n 



Figure 1 . A schematic representation of the evolution of singlet number density and how it is divided 
between the condensate and excitations. The dashed line shows the constant comoving number 
density, the purple curve shows the number density of singlet particles and the blue one shows the 
number density of the condensate. The figure corresponds to the case where the interactions are 
sufficiently strong to allow for the condensate to decay away (T ~ H) before the quadratic part of the 
singlet potential starts to dominate the evolution (at T trans ). If the condensate survives until T trans , 
it practically does not decay at all and its comoving number density remains almost constant. 


effectively relativistic matter until /c*/a ~ m s and they become non-relativistic. If the decay 
of the condensate is complete, the corresponding transition temperature T tra ns is given by 
Eq. (3.3) as before. For T < T trans the energy density of singlet particles is thus given by 
p s = m s n s which yields 



for the present abundance of singlet dark matter. 

The parametric dependence of the dark matter yield is illustrated in Figures 2 and 3. 
The red domain in the figures marks the regime where the condensate decays completely 
and dark matter consists of singlet particles. This domain exists only for sufficiently large 
values of the singlet self-coupling A s or the portal coupling A s h- For smaller couplings the 
condensate decay never becomes efficient and the final dark matter component will consist 
of the singlet condensate instead of particles. 

As reviewed in Sec. 2, for superweak portal coupling the correct dark matter abundance 
can be generated at low temperature [27]. Here we have demonstrated that the primordial 
singlet condensate crucially affects the abundance and properties of singlet dark matter 
providing non-trivial boundary condition for the low temperature dark matter production. 
We leave for future work the careful matching of the low temperature particle production 
with the high temperature freeze-in we have described here. 

5.2 Case II: The condensate survives 

After the transition from quartic to quadratic oscillations the decay rate of the singlet con¬ 
densate decreases faster than the Hubble rate. In the parameter range where r so _> ss <C H 
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until the end of the quartic epoch the condensate therefore practically does not decay at all 
but is merely redshifted according to Eq. (3.3). 

In more detail, the corresponding Boltzmann equation for singlet particles (4.5) in the 
quartic regime T > T trans is given by, 


,;y( 4 ) JV I'TOhN p( 4 ) p(4) 

_ A 1 l T ) L h-*sh _ 1 s o^ ss 

dT K 2 (^) Hs h T h Hs h T [ 

where Y s = n s /sb denotes the singlet number density normalized by the entropy density of 
the bath Sb and where we used T ~ — HT , which is an excellent approximation above the 
EW scale. With the rates given in Eqs. (4.7), (4.9) and (4.10), the solution of Eq. (5.8) is 



After the transition to the quadratic regime T < Tt rans the kinematical suppression 
renders the singlet particle production negligible and the comoving particle number freezes 
to a constant value. The corresponding present particle abundance is 

(§§) = 2-286 x 10 9 (^) T s ( 4 ) (To) • (5.10) 

The final yield corresponding to the correct DM abundance today is depicted in Figures 
2 and 3 for representative values of model parameters. The regions below the red domains 
correspond to the case where the condensate survives. The purple areas just over the blue 
regions depict the case II, where there is a small fraction of relic density also in the singlet 
particles. The dotted line denotes the upper bound for an isocurvature DM component, 
allowing only ~1% of non-thermally generated DM contributing to the observed total DM 
abundance [9]. The dashed lines, from thinnest to thickest denote 5%, 20%, 50% and 80% 
abundances, respectively. For example, in Figure 2 the upper boundary of the red region 
corresponds to 0.1% abundance and the lower boundary of the blue band to 100%. The 
yellow vertical band shows the standard freeze-in scenario [27, 34, 35], where only the low 
temperature processes (h ss ) produce 1-100% of the present DM abundance, assuming 
f s = 0 at T = Tew- Note that the portal coupling is always required to be A s h < 10 ' in 
order to avoid singlet thermalization before the electroweak scale. 

As we have discussed above, one important consequence of the primordial singlet con¬ 
densate is that it will lead to isocurvature perturbations which are heavily constrained by 
their imprints on the CMB anisotropies [9]. For the specific SM extension (2.1) we have con¬ 
centrated on here, this constrains the singlet dark matter component to constitute at most 
only a very small fraction of the total dark matter abundance, as shown in Figures 2 and 3. 

However, extending beyond the simplistic scenario (2.1) these constraints could be eas¬ 
ily relaxed. In particular, we can entertain a thought that there are additional fields and 
interactions in the portal sector. These should not affect the high temperature evolution we 
have considered here, but could provide additional channels for the singlet sector to inter¬ 
act with the SM fields at temperatures below the electroweak transition. These interactions 
could convert the isocurvature modes into the observed adiabatic perturbations. A concrete 
possibility would be a superheavy scalar attaining a non-zero vev only at low temperatures. 
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Figure 2. The final dark matter yield, shown in terms of the portal coupling A s h and the self¬ 
interaction coupling A s . Different regions where the model components (singlet condensate, singlet 
particles) produce 0.1-100% of today’s dark matter abundance are shown in red (particles only), 
blue (condensate only) and purple (both particles and condensate); see the main text for further 
explanations. The slashed contours refer to 80%, 50%, 20% and 5% abundances, listed from thickest 
to thinnest contour. The green region shows where the singlet would constitute more than 100% of 
today’s DM abundance. The dotted line denotes the upper bound (~1%) for an isocurvature DM 
component. The vertical yellow band shows the standard scenario where only the low temperature 
processes produce 1-100% of today’s DM abundance. 

6 Conclusions and outlook 

We have considered in detail the consequences of inflationary initial conditions on the dynam¬ 
ics of (extended) scalar sectors with light excitations. Concretely, we considered the simple 
benchmark model where the dark matter is constituted by a Z 2 symmetric real scalar field. 
Our main results are equations (5.7) and (5.10) together with Figures 2 and 3. They have 
strong implications on presently popular models of dark matter production via the freeze-in 
mechanism. 

One particularly important feature we have uncovered is that, contrary to what is often 
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Figure 3. The final dark matter yield, shown in terms of the singlet mass m s and tensor-to-scalar 
ratio r (left panel) or the self-interaction coupling A s (right panel). The different colours and contours 
are the same as in Figure 2. 


assumed, A s is not irrelevant for the production of the dark matter abundance. Another im¬ 
portant feature is that significant particle production via the freeze-in mechanism is possible 
at high-temperatures, above the electroweak scale, due to primordial condensates following 
from the initial conditions set by inflation. Consequently, the initial conditions for the singlet 
abundance to be negligible at T ~ Tew ma Y not be valid. More generally, our analysis has 
revealed a novel connection how fundamental inflationary physics is imprinted on the dark 
matter abundance. 3 

Within the singlet scalar extension of SM we have shown how the freeze-in scenario 
severely constrains both the self-interaction coupling A s and portal coupling A s h- Our results 
can be also extended to provide new constraints on models in which the scalar acts only as 
a mediator and decays further to the actual DM particle such as a sterile neutrino. We will 
consider scenarios beyond the simple benchmark one treated in this paper in a forthcoming 
work. 
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